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experiencing  low  attenuation.  They  have  been  successfully  used  for  damage  detection  in  structures  of  relatively  low 
geometric  complexity  such  as  plates  and  cylindrical  pipes.  The  use  of  guided  waves  for  this  purpose  becomes 
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ABSTRACT 

Structural  health  monitoring  (SHM)  technologies,  which  use  integrated  sensing  for  damage  detection,  are  expected  to 
improve  system  reliability,  availability,  and  operational  cost.  Guided  waves  can  propagate  great  distances  while 
experiencing  low  attenuation.  They  have  been  successfully  used  for  damage  detection  in  structures  of  relatively  low 
geometric  complexity  such  as  plates  and  cylindrical  pipes.  The  use  of  guided  waves  for  this  purpose  becomes 
increasingly  difficult  as  the  geometric  complexity  of  the  structure  increases.  Aerospace  structural  components  such  as 
fuel  tanks,  wings,  etc.  often  are  comprised  of  substructures  that  consist  of  plates  with  integral  stiffeners.  This  work 
reports  on  finite  element  simulations  of  guided  waves  in  integrally  stiffened  plate  structures.  In  these  studies,  the  guided 
waves  are  generated  by  PZT  wafer-type  transducers  mounted  on  the  structure.  Transient  dynamic  finite  element 
simulations  using  PZFlex,  in  2D  and  in  3D,  were  used  to  model  both  the  structure  and  transducers.  The  interaction  of  the 
guided  waves  with  cracks,  simulated  by  notches  of  varying  dimensions,  is  also  modeled.  This  allows  appraisal  of  the 
sensitivity  of  various  modes  for  crack  detection  by  providing  insight  into  mode  conversion  and  scattering  resulting  from 
the  guided  wave  and  crack  interaction. 
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Flealth  Monitoring 


1.  INTRODUCTION 


Structural  Flealth  Monitoring  (SF1M)  has  the  functional  goals  of  determining  when,  where,  and  to  what  extent  structural 
damage  has  occurred  through  the  use  of  integrated  sensing.  A  possible  further  functional  objective  is  the  sensing  of 
precursors  to  damage.  The  driving  operational  forces  behind  the  developing  field  of  Structural  Flealth  Monitoring 
(SF1M)  are  the  goals  of  increasing  system  reliability,  improving  system  availability,  and  reducing  inspection  and 
maintenance  costs.  The  use  of  integrated  ultrasonic  transducers  and  guided  waves  is  a  technology  that  is  viewed  as 
having  great  potential  for  SF1M. 

Guided  waves  have  the  ability  to  cover  a  large  region  of  a  structure  from  a  transducer  or  group  of  transducers  resident  in 
a  small  area.  Much  of  the  research  and  implementation  performed  to  date  has  focused  on  inspection  of  plate  and  tubular 
structures[l].  This  is  reasonable  in  view  of  the  relatively  low  geometric  complexity  of  these  structures.  Another  area  that 
has  received  considerable  attention  is  multilayer  structures.  Lap  joints  and  ‘doublers’  are  examples  of  typical  aerospace 
applications^]. 

An  integrally  stiffened  structure  results  when  a  plate  has  a  stiffener  attached,  usually  at  a  right  angle  to  the  plate  surfaces. 
The  simplest  example  of  this  geometry  is  the  extruded  T-Section.  The  cross  section  of  this  geometry  is  illustrated  in 
Figure  1.  They  represent  a  somewhat  greater  geometric  complexity  than  plate,  tube,  and  multilayer  structures.  Integrally 
stiffened  structures  have  many  applications  in  aerospace  vehicles.  They  commonly  appear  as  substructures  in  wings, 
airframes,  and  fuel  tanks. 
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Coupling  of  Lamb  modes  at  a  T  junction  was  experimentally  investigated  by  Morvan,  et  al.  [3].  Greve,  et  al.  reported  on 
a  numerical  simulation  of  Lamb  wave  propagation  in  a  plate  girder(“I-Beam”)  geometry[4].  This  simulation  had  an  So 
Lamb  wave  generated  in  the  web  and  studied  the  interaction  of  the  S0  mode  with  a  crack  in  the  weld  bead  on  the  web 
side  of  the  web-flange  interface.  A  two-dimensional  T-Section  geometry  with  fillets  with  and  without  notch  defects  is 
numerically  simulated  by  Aldrin,  et  al.[5]. 

The  objective  of  this  study  is  to  investigate  the  generation  and  interaction  of  guided  waves  in  a  three  dimensional  T- 
section  geometry  both  with  and  without  notches  to  simulate  cracks.  The  waves  are  generated  in  the  base  of  the  T- 
Section,  which  would  correspond  to  a  girder  flange.  This  modeling  includes  an  explicit  model  of  a  bonded  transducer 
disk  to  more  closely  simulate  experimental  observation  and  also  allow  observation  of  the  three  dimensional  effects. 

2.  FINITE  ELEMENT  MODEL 

In  this  work  we  examine  the  two  dimensional  and  three  dimensional  T-Section  problems.  The  2D  geometry 
corresponding  to  the  3D  cross-sectional  geometry  is  shown  in  Figure  1.  The  model  is  geometrically  parameterized  in 
terms  of  the  thickness  and  length  of  the  left  leg,  thickness  and  length  of  the  vertical  leg,  thickness  and  length  of  the  right 
leg,  and  the  extrusion  width.  The  bonded  transducer  is  also  parameterized  in  terms  of  the  transducer  radius  and 
thickness.  The  transducer  location  is  parameterized  in  terms  of  offsets  from  the  symmetry  planes  of  the  3D  T-Section. 
There  are  also  validity  checks  which  keep  the  sensor  inside  the  bounds  of  the  part  and  do  not  allow  placement  too  close 
to  the  vertical  leg.  The  parameter  values  of  the  model  for  this  study  are:  f  =  2.21  mm,  t2  =2.21  mm,  t3  =2.21  mm,  tt  = 
0.19  mm,  Wi  =30.0  mm,  w2  =  90.0  mm,  w3  =  12.7  mm,  wt  =80.0  mm,  and  rt  =  5.0  mm. 

The  3D  geometry  of  the  model  showing  the  transducer  and  the  length,  d  =  30.0  mm,  in  the  extrusion  direction  is  shown 
in  Figure  2.  The  transducer  disk  bonded  to  the  top  side,  centered  on  the  extrusion  dimension  symmetry  plane,  of  the 
plate  is  made  of  PZT-5H.  The  transducer  electrodes  are  on  its  top  surface  and  bottom  surface.  The  signal  is  applied  at  the 
top  surface  with  the  bottom  surface  being  grounded.  The  interface  between  the  transducer  and  the  structure  is  modeled  as 
being  ‘welded’,  i.e.  the  adhesive  layer  is  not  explicitly  modeled  and  the  displacement  and  force  are  continuous  across  the 
interface.  The  piezoelectric  equations  are  solved  within  a  volume  of  a  rectangular  parallelepiped  (‘electric  window’) 
whose  dimensions  are  equal  to  the  diameter  and  thickness  of  the  transducer.  The  crack  location  is  also  shown  in  Figure 
1.  The  notches  at  the  base  were  all  0.29mm  wide(longitudinal  direction).  Three  notch  depths(transverse  direction)  were 
modeled  in  addition  to  the  notch  free  configuration.  The  notch  depths,  tn,  were  0.402,  0.804,  and  1.206  mm.  The  notches 
in  the  3D  model  ran  across  the  entire  extrusion  depth,  30  mm. 

The  voltage  applied  to  the  transducer  in  the  model  uses  a  single  cycle  sine  function  at  a  frequency  of  500MHz.  The 
signal  is  10  volts  peak  to  peak.  The  outer  planes  have  an  absorbing  boundary  condition  applied.  This  removes  the 
complications  of  reflected  signals  and  the  model  can  be  viewed  as  a  piece  of  a  much  larger  structure.  The  f  d  product  is 
1.11  mm  MHz  and  is  in  the  region  of  the  Lamb  wave  dispersion  curve  where  only  the  A0  and  So  modes  are  expected  to 
be  present[6-7].  The  mesh  density  in  the  model  was  approximately  60  elements  per  wavelength  at  500  MHz.  This  model 
was  generated  and  run  using  the  PZFlex  Finite  Element  software[8]. 
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Fig.  2.  Three  Dimensional  T-section  (Integrally  Stiffened  Plate)  showing  extrusion  depth  and  transducer  and  location. 


3.  DISCUSSION  OF  RESULTS 


The  response  of  the  2D  and  3D  T-Sections  to  the  transducer  input,  both  without  and  with  notches,  were  run  for  this 
study.  The  stresses  displayed  in  the  figures  are  the  difference  between  the  normal  stresses  in  the  directions  parallel,  x, 
and  orthogonal,  y,  to  the  centerline  of  the  bottom  plate.  This  stress  is  what  is  observed  in  photoelastic  wave  visualization 
experiments[9].  The  terminology  of  ‘photoelastic’  stress  is  adopted  for  this  paper.  This  can  be  expressed  as: 


<Jp  Gyy  -  Gxx  ( 1 ) 

where  ap  is  the  photoelastic  stress,  ayy  is  the  y  normal  stress,  and  crvv  is  the  x  normal  stress.  These  stresses  are  calculated 
through  the  use  of  the  scripting  environment  available  in  PZFlex  which  allows  access  to  the  data  structures  which 
contain  the  native  field  variables  that  the  software  does  calculate. 
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3.1  2D  Model 


The  photoelastic  stress,  for  the  un-notched  2D  model,  is  shown  in  Figure  3  for  a  time  shortly  after  the  wave  front  of  the 
response  has  encountered  the  T  junction.  This  figure  contains  vertical  bands  on  the  left  hand  side  which  are  indicative  of 
the  S0  mode.  The  So  mode  has  a  higher  propagation  velocity  than  the  A0  mode  and  is  separating  out  of  the  mixture  of  the 
two  modes  generated  by  the  transducer  pulse.  The  patterns  closer  to  the  transducer  have  vertical  bands  which  change 
sign  but  at  a  point  away  from  the  centerline.  This  indicates  a  mixture  of  both  A0  and  S0  modes.  The  bottom  view  in  the 
figure  is  a  close-up  of  the  T  junction.  The  interesting  phenomenon  observed  in  this  view  is  that  the  pattern  induced  in  the 
vertical  stiffener  appears  to  be  indicative  of  an  A0  mode  even  though  only  excitation  is  the  ‘pure’  S0  wave  front.  A 
possible  explanation  for  this  is  that  the  SO  wave,  of  frequency  f,  passing  by  the  stiffener  excites  it  at  that  frequency  f. 
Since  the  wave  is  approximately  sinusoidal  in  the  thickness  direction  the  excitation  will  consist  of  spatially  symmetric 
and  anti-symmetric  portions.  The  stiffener  thickness  is  the  same  as  base  plate  thickness  which  results  in  the  same  fd 
value  as  the  base  plate. 
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Fig.  3.  Photoelastic  stress,  Pa,  response  of  a  2D  T-section  to  a  single  cycle  sine  pulse  transducer  input  at  time  =  14.94  ps. 
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Fig.  4.  y  displacement  profiles  at  various  times  along  top  of  plate  for  a  2D  T-section  subjected  to  a  single  cycle  sine  pulse 
transducer  input. 


The  transverse,  y,  displacements  along  the  top  of  the  plate  on  the  center  cross-section  at  various  times  are  shown  in 
Figure  4.  These  show  the  S0  modes  separating  from  the  A0  modes  as  time  progresses.  The  plot  for  14.9  ps  (1.49  e-05 
seconds)  corresponds  to  the  photoelastic  stress  plot  in  Figure  3.  The  differences  in  the  baseline  (i.e.  un-notched)  and 
notched  cases  are  difficult  to  differentiate  visually.  Figure  5  is  a  plot  of  the  y  displacement  at  various  x  locations  on  top 
of  the  plate.  The  propagation  of  the  wave  front  through  the  structure  may  be  clearly  seen.  Figure  6  is  a  plot  of  the 
differential  signal  between  the  baseline  and  notched  y  displacement  at  various  x  locations  on  top  of  the  plate.  The 
differential  signal  is  observed  to  originate  at  the  notch  location  and  propagate  out  in  both  directions.  The  differential 
signals  of  the  three  different  notch  sizes  received  at  the  x  location  value  of  -235  mm  are  shown  in  Figure  7.  These  signals 
appear  to  increase  in  magnitude  and  experience  an  increasingly  lagging  phase  shift  with  increasing  notch  depth.  This 
could  be  the  result  the  reduction  in  plate  thickness  causing  a  shift  along  the  dispersion  curve  in  the  notch  region.  Further 
investigation  of  this  is  desirable. 
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Fig.  7.  Differential  y  displacement  of  0.4  mm,  0.8  mm,  1.2  mm  and  base  cases  versus  time  on  top  of  plate  at  x  location  -235 
mm  for  a  2D  T-section  subjected  to  a  single  cycle  sine  pulse  transducer  input. 

3.2  3D  Model 

Figure  8  shows  the  photoelastic  stresses  along  cross  sections  of  the  3D  geometry.  The  cross  sections  are  taken  normal  to 
the  ‘extrusion’  direction  of  the  T-section.  There  are  five  cross  sections,  taken  in  equal  increments,  starting  at  the  center 
of  the  T-section  and  ending  at  the  outer  edge.  Since  the  model  is  symmetric,  only  half  of  the  model  is  shown.  The  sub¬ 
windows  show  close-ups  of  the  T  junction  where  the  vertical  stiffener  meets  the  base  plate.  These  start  at  the  center  and 
move  outward  in  equal  increments  toward  the  edge.  A  close-up  of  the  outer  boundary  is  not  shown. 

The  character  of  the  response,  especially  along  the  center  cross  section  of  the  model,  displays  great  similarity  to  the  2D 
case.  The  transducer  in  the  3D  model  is  a  circular  disk.  The  axisymmetry  of  the  transducer  geometry  gives  the 
expectation  that  circular  crested  Lamb  waves[10]  will  be  generated  instead  of  the  straight  crested  Lamb  waves  that  result 
in  the  two  dimensional  plane  strain  case.  Circular  crested  waves  have  displacement  functions  through  the  thickness  that 
are  similar  to  the  straight  crested  case  but  they  are  weighted  by  Bessel  functions  in  the  radial  direction.  The  wave  lengths 
are  also  functions  of  the  radial  coordinate  but  asymptotically  approach  the  straight  crested  case.  In  addition,  the  straight 
crested  dispersion  curves  also  display  great  similarity  to  the  circular  crested  case.  The  subsequent  discussion  of  S0  and 
A0  modes  refers  to  the  circular  crested  equivalents  of  the  straight  crested  Lamb  modes.  Figure  8  indicates  a  mixture  of 
A0  and  So  modes  analogous  to  the  2D  situation.  The  close-up  of  the  3D  T  junction  also  displays  an  A0  type  pattern 
induced  in  the  vertical  stiffener  even  though  it  has  only  experienced  excitation  by  a  ‘pure’  So  wave  front. 

Figure  9  displays  the  transverse,  y,  displacements  along  the  top  of  the  plate  on  the  center  cross-section  at  various  times. 
The  S0  modes  appear  to  separate  from  the  A0  modes  as  time  progresses.  The  plot  for  1 .49  e-05  seconds  is  close  to  the 
1.64e-05  seconds  of  the  photoelastic  stress  plot  in  Figure  8.  A  plot  of  the  y  displacement  at  various  x  locations  on  top  of 
the  plate  is  shown  in  Figure  10.  The  wave  front  propagation  through  the  structure  is  clearly  shown.  The  difference 
between  the  un-notched  (base)  and  the  notched  cases  is  hard  to  detect  visually.  A  plot  of  the  differential  signal  between 
the  baseline  and  notched  y  displacement  at  various  x  locations  on  top  of  the  plate  is  illustrated  in  Figure  1 1.  As  in  the  2D 
situation,  the  differential  signal  is  observed  to  originate  at  the  notch  location  and  propagate  out  in  both  directions.  The 
differential  signals  resulting  from  the  three  different  notch  sizes  received  at  the  x  location  value  of  -235  mm  are  shown 
in  Figure  12.  These  signals  appear  to  increase  in  magnitude  and  experience  an  increasing  phase  shift  with  increasing 
notch  depth  in  a  manner  analogous  to  that  observed  in  2D.  Comparison  of  Figures  4  and  9  shows  the  greater  attenuation 
of  the  3D  guided  waves.  This  attenuation  is  also  evidenced  by  the  greater  magnitude  of  the  2D  differential  signal  in 
Figure  7  versus  the  3D  signal  in  Figure  12. 
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Response  of  an  T- Section  to  a  1  cycle  sin  pulse  transducer  input  absorbing  BCs  0.5MHz 


Fig.  8.  Photoelastic  stress.  Pa,  response  of  a  3D  T-section  to  a  single  cycle  sine  pulse  transducer  input  at  time  =  16.23  ps. 


Fig.  11.  Differential  y  displacement  comparing  base  and  1.2  mm  cases  versus  time  on  top  of  plate  at  various  x  locations  for 
a  3D  T-section  subjected  to  a  single  cycle  sine  pulse  transducer  input. 
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Fig.  12.  Differential  y  displacement  of  0.4  mm,  0.8  mm,  1.2  mm  and  base  cases  versus  time  on  top  of  plate  at  x  location  - 
235  mm  for  a  3D  T-section  subjected  to  a  single  cycle  sine  pulse  transducer  input. 


4.  CONCLUSIONS 


The  response  of  a  2D  and  3D  integrally  stiffened  geometry  with  and  without  notching  to  excitation  by  a  piezoelectric 
transducer  was  modeled  by  using  Finite  Element  analysis.  The  response  along  the  center  section  of  the  3D  structure  to  a 
circular  disk  transducer  is  quite  similar  in  character  to  the  2D  model  excited  by  a  strip  transducer.  In  comparison  to  the 
straight  crested  Lamb  waves  generated  in  the  2D  case,  the  circular  crested  Lamb  waves  of  the  3D  case  do  attenuate  in 
perfectly  elastic  material  due  to  spreading.  This  does  not  present  a  problem  if  the  distances  involved  are  not  great  but 
large  enough  that  the  waves  have  asymptotically  approached  straight  crested  behavior.  Both  models  had  the  appropriate 
SO  and  A0  modes  generated.  Separation  between  the  two  modal  types  was  observed  as  propagation  progressed  due  to 
their  different  velocities. 

The  change  in  response  to  the  excitation  due  to  different  sized  notches  at  the  base  of  the  T  junction  was  also  modeled. 
The  comparison  is  done  with  the  baseline  (un-notched)  case.  There  is  a  differential  observed  in  the  displacements  along 
the  centerline  of  the  top  of  the  plate  in  the  signals  received  with  and  without  the  notch  present.  These  differential  signals 
display  magnitudes  and  phase  shifts  that  scale  with  increasing  notch  size. 

Several  possibilities  present  themselves  as  areas  for  future  work.  A  comparison  with  experimental  observation  is 
desirable.  A  study  of  the  case  where  a  receiving  transducer  is  explicitly  modeled  and  the  size  of  notch  required  for 
detection  would  further  tie  this  work  to  actual  practice.  Another  valuable  investigation  is  the  effect  of  the  stiffener 
geometry  on  the  mode  conversion  and  the  effect  of  reflection  from  the  top  surface  of  the  stiffener  back  into  the  base 
plate. 
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